Monte Carlo simulations of the structure of Pt-based bimetallic nanoparticles 
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Pt-based bimetallic nanoparticles have attracted significant attention as a promising replacement 
for expensive Pt nanoparticles. In the systematic design of bimetallic nanoparticles, it is impor- 
tant to understand their preferred atomic structures. However, compared with unary systems, alloy 
nanoparticles present more structural complexity with various compositional configurations, such 
as mixed-alloy, core-shell, and multishell structures. In this paper, we developed a unified empirical 
potential model for various Pt-based binary alloys, such as Pd-Pt, Cu-Pt, Au-Pt, and Ag-Pt. Within 
this framework, we performed a series of Monte Carlo (MC) simulations that quantify the energet- 
ically favorable atomic arrangements of Pt-based alloy nanoparticles: an intermetallic compound 
structure for the Pd-Pt alloy, an onion-like multi-shell structure for the Cu-Pt alloy, and core-shell 
structures (Au@Pt and Ag@Pt) for the Au-Pt and Ag-Pt alloys. The equilibrium nanoparticle 
structures for the four alloy types were compared with each other, and the structural features can 
be interpreted by the interplay of their material properties, such as the surface energy and heat of 
formation. 

PACS numbers: 61.46.+W, 36.40.Ei, 64.70.Nd 



o 

>\ 

<^ 

(N 
O 



X: 



INTRODUCTION 



Due to their high surface-to-volume ratio or quan- 
tum effects at extremely small sizes, metal nanoparti- 
cles exhibit unique chemical and physical properties that 
are distinct from bulk materials''-. Among the many 
applications of nanoparticles, such as electronic, opti- 
cal, and magnetic devices^r"^, their electro-catalytic ef- 
fects are the most beneficial, especially for many effec- 
tive energy conversion techniques, e.g., the performance 
of proton-exchange membrane fuel cells (PEMFCs) and 
direct methanol fuel cells (DMFCs) as power resources 
for electric devices^"—. 

Platinum (Pt) is practically the only catalyst active for 
hydrogen oxidation, methanol oxidation and oxygen re- 
duction at low temperature and it is one of the most effi- 
cient metal catalysts for fuel cells and similar energy con- 
version techniques^. As mentioned in many articles^^— , 
Pt nanoparticles with sizes of a few nanometers provide 
exclusively outstanding catalytic activity/stability. The 
primary disadvantage of using Pt nanoparticles in appli- 
cations is the high cost of Pt (which is one of the main 
obstacles to the commercialization of fuel cells). Two so- 
lutions have been suggested for reducing the industrial 
demand for Pt in fuel cells: replacing Pt with alter- 
native non-noble catalysts and reducing the Pt-loading 
by exploiting non-precious supports^. In this regard, 
the design of Pt-based bimetallic nanoparticles to re- 
place Pt nanoparticles has attracted significant atten- 
tion^riS. For the systematic design of bimetallic nanopar- 
ticle catalysts, it is very important to understand their 
stable atomic structures because the chemical and phys- 
ical properties of various Pt-based alloy metal nanopar- 
ticles are highly dependent on their size, shape, compo- 



sition, morphology, and structural stability. 

Compared with unary catalysts, alloy catalysts 
present more structural complexity because the two 
components can exhibit various structural configura- 
tions/modifications. For example, bimetallic nanoparti- 
cles can exhibit the structures of ordered/random mixed- 
alloj«^, core-shelliiii^, and multishell nanoparticlesi^ii^. 
The equilibrium structure of bimetallic nanoparticles is 
also sensitive to the alloy components, involving the in- 
terplay of several material properties, such as surface en- 
ergy, heat of mixing, and strain effect due to the atomic 
size difference. Because of these complex factors, it is 
difficult to predict the most stable structure of bimetal- 
lic nanoparticles and their thermodynamic properties. 
In this paper, we developed a unified embedded atom 
method (EAM) potential model for simulating the struc- 
tures of various Pt-based binary alloys, such as Pd-Pt, 
Cu-Pt, Au-Pt, and Ag-Pt. Within the framework of our 
potential model, we performed a series of Monte Carlo 
(MC) simulations to quantify the energetically favorable 
atomic structures of Pt-based alloy nanoparticles. 



II. INTERATOMIC POTENTIAL MODEL 

A. Embedded atom method potentials for Pt, Pd, 
and novel metals (Ag, Au, Cu) 

The core ingredient of atomic-scale simulations is how 
to evaluate the potential energy of systems and the in- 
teratomic forces as a function of the positions of the 
atoms. Although it originates from the electronic struc- 
tures, in classical atomistic simulations, such as molec- 
ular dynamics (MD) and Monte Carlo (MC) methods. 
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the interaction between atoms is usually described by 
empirical potentials, and therefore the reliability of the 
simulation result is entirely dependent on the reality and 
accuracy of the interatomic potentials. For most metallic 
systems, the embedded atom method (EAM) potential 
is widely used, and several potentials have been devel- 
oped for elemental Pt and some other face-centered cubic 
(fee) alloys^^"^^. Especially in many applications of al- 
loy systems, EAM potential models of analytic functional 
form have been frequently used because of their simplic- 
ity and extendibility to multi-component systems. Here, 
we adopted EAM potential models for fee metals that 
were originally developed by A. Voter— and later modi- 
fied by us^S. The total energy of the system is given by 
the usual EAM form: 



(1) 



where Fs- {p) is the energy associated with embedding 
atom of type Si in a uniform electron gas of density p 
and </>£,. _ sj (r) is a pairwise interaction between atoms of 
type Si and Sj separated by a distance r. The electron 
density is given by 



Pi 



(2) 



where the atomic electron density function /(r) is taken 
as the density of a hydrogenic 4s orbital: 



(3) 



Here, /3 is an adjustable fitting parameter that quantifies 
the distance over which the electron density decays away 
from an atom position and /o is a prefactor. (It can be an 
arbitrary value for a unary system because it cancels out 
when combined with the embedding function, but it is an 
adjustable parameter that should be determined for an 
alloy system.) We chose /o to be l/{J{reqNist) for con- 
venience, so that the electron density at the equilibrium 
crystals is approximately normalized to unit value. 

Within the formalism of A. Voter—, the pair potential 
term, (j>{r) , is chosen to take a Morse potential form with 
a minor additional term: 
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where Dm, Rm, and um are adjustable fitting pa- 
rameters that define the well-depth of the Morse pair- 
interaction, the position of the minimum, and a measure 
of the curvature at the minimum, respectively. The sec- 
ond term was introduced to tune the melting point of the 
EAM model: the parameter 6 represents the magnitude 
of the pair-interaction tuning between the second and 



third nearest neighbor positions, r2 and r^, with the aid 
of the Heaviside step- function, 9{r). Details regarding 
the meaning of this term will be provided elsewhere^S. 

Finally, the embedding function was numerically de- 
termined so that the total energy of the reference system 
as a function of dilation satisfies the following universal 
binding energy relatioiii^i: 
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(5) 

with a — y^QBn/ Eq, where a is the dilated lattice con- 
stant, ao is the equilibrium lattice constant and Eq, B, 
and 57 are the cohesive energy, bulk modulus, and equi- 
librium atomic volume of the reference lattice, respec- 
tively. The potential interactions were smoothly cut 
off at r = Vcut (usually between the third- and fourth- 
nearest-neighbor shells of a static fee crystal) to ensure 
that the interatomic potential and its first derivatives are 
continuous. 

Because it takes such an analytic functional form, this 
potential automatically gives exact matching with the 
experimental values of the equilibrium lattice parameter 
oq, cohesive energy Eq, and bulk modulus B via Eq. ([S]). 
Thus, we need to fit the remaining six adjustable param- 
eters {D]^j , RmjChm, and r^ut) for a single compo- 
nent, and those are obtained by optimization relative to 
some material properties. The parameters in the Voter's 
potentials were fitted to a series of basic properties, in- 
cluding the three cubic elastic constants (Cn, C12, and 
C44) and the vacancy formation energy, E^aci of the face- 
centered cubic (fee) crystal for Pt, Pd, and the novel 
metals^S,. Here, instead of using the original parameter 
set suggested by Voter—, we re-optimized the potential 
parameters (except Vcut, which is of less importance), by 
including the surface energy in the target properties^. 
The resulting parameter sets are shown in Table |T] and 
the calculated material properties are compared with the 
target values in Table HIl 

Among the various materials properties, we believe 
that the surface energies of the elements are of partic- 
ular importance for the simulation of nanoparticles. Fig- 
ure [T] shows the surface energies of the five elements (Pt, 
Pd, Cu, Au, and Ag) that were calculated by our EAM 
potential compared with first-principle calculations^^"— 
and experimental data3ii^. While there is a diverse data 
distribution in the first-principle calculation results, the 
order of the experimental surface energy was consistently 
7Pt > iPd > iCu > lAu > lAg- The surfacc ener- 
gies of the atomic species are expected to be one of the 
most influential properties in the energetics of nanopar- 
ticles due to the high surface-to-volume ratio, and we 
believe that the relative value of the surface energy will 
play a role in determining the most favorable nanopar- 
ticle structures and their structural behavior. Note that 
the present EAM potential model predicts the correct 
order of the surface energies among the elements even 
better than the first-principle calculations and hence is 
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TABLE I: Fitted potential parameters for five fee metals: Pt, 
Pd, Cu, Au, and Ag. 



Parameter 


Pt 


Pd 


Cu 


Au 


Ag 


p (A-) 


3.580 


3.400 


3.990 


3.795 


3.960 


au (A-i) 


1.779 


1.501 


1.838 


1.753 


1.685 


Dm (eV) 


0.783 


1.682 


0.781 


0.695 


0.693 


Rm (A) 


2.577 


2.343 


2.321 


2.537 


2.538 


rcut (A) 


5.576 


5.412 


4.961 


5.516 


5.542 


5 (eV) 


0.01 


0.06 


0.0 


0.049 


0.0235 



>2 2000 



§ 1500 



3 
t/J 



■ Y^^^^_^^^, (Experiment) 
■I'loo I'lio' ■*■ ■■'ill (PresentEAM) 
^100 ■'no ^ ^111 (DFT calculations) 
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TABLE IL Materials properties used in fitting for fee Pt, Pd, 
Cu, Au, and Ag. "Where two numbers are given, the top 
number is the calculated value predicted by the present EAM 
model, and the lower number in parentheses is the experimen- 
tal value from Ref.[23 (elastic properties) or the first-principle 
calculation results from Refs. 23-26 (surface energies). 



Properties 


Pt 


Pd 


Cu 


Au 


Ag 


ao (A) 


3.924 


3.891 


3.615 


4.078 


4.085 


Ecoh (e'V/atom) 


5.77 


3.91 


3.54 


3.93 


2.85 


B (GPa) 


283 


195 


142 


167 


104 


Cii (GPa) 


319 


234 


180 


185 


121 




(347) 


(234) 


(176) 


(186) 


(124) 


Ci2 (GPa) 


264 


175 


124 


158 


94.9 




(251) 


(176) 


(125) 


(157) 


(93.4) 


C44 (GPa) 


77 


70.5 


79.8 


40.3 


43.7 




(77) 


(71.2) 


(81.8) 


(42.0) 


(46.1) 


(eV) 


1.4 


1.58 


1.3 


1.0 


1.2 




(1.5) 


(1.54) 


(1.3) 


(0.9) 


(1.1) 


7100 (mJ/m^) 


1891 


1815 


1435 


1038 


1009 




(1889) 


(1900) 


(1450) 


(1444) 


(1200) 


7110 (mJ/m^) 


2079 


2012 


1605 


1138 


1120 




(1920) 




(1530) 


(1700) 


(1290) 


7111 (mJ/m^) 


1580 


1602 


1266 


865 


886 




(1531) 


(1880) 


(1300) 


(1040) 


(1120) 



expected to be transferable to nanoparticle systems with 
a large surface-to- volume ratio. 

Figure [2] shovirs the melting temperature of the five el- 
ements predicted by our EAM potential model. (The 
melting points of the fee metals 'were obtained via mi- 
crocanonical ensemble molecular dynamics simulations of 
the solid- liquid coexistenceSS.) With the aid of the tun- 
ing term in the pair-potential function, our EAM model 
reproduces the experimental melting temperature ■within 
less than a 5% difference. 



B. Interatomic potentials for Pt-based alloys 

Significant effort has also been made to develop empir- 
ical potential models for metallic alloys. For alloy poten- 



Pt Pd Cu Au Ag 

Element 

FIG. 1: Surface energies of the five metals: Pt, Pd, Cu, Au, 
and Ag. The calculated values based on our EAM model are 
compared with the first-principle calculations and the experi- 
mental data. The experimental surface energies (solid circles) 
correspond to averaged values over the various surfaces, which 
are from Refs. [23 and [2^. The calculated values are mea- 
sured on the (100), (110), and (111) surfaces (Refs. [US). 
The first-principle calculations show diverse data depending 
on their exchange-correlation functionals. 
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FIG. 2: Melting temperatures of the five metals (Pt, Pd, Cu, 
Au, and Ag): Our EAM model reproduces the experimental 
melting temperatures to within less than 5%. 



tials, ■we describe the interaction bet^ween atoms of differ- 
ent species, which are often called the cross-interaction 
potentials. In the early versions of EAM potentials tak- 
ing analytic forms, cross-interaction potentials were pro- 
posed to take a simplified combination of the individual 
pair-interactions based on the potentials for pure met- 
als, e.g., the cross-interaction potential was assumed to 
be the geometric mean of the monatomic pair poten- 
tialsi^ or to be a density-weighted combination of the 
monatomic pair potentialaiii^. This approach was fairly 
reasonable for several binary solid solutions of fee met- 
als, but not sufficient for a quantitative study of the alloy 
properties. Currently, even for a system in which no ad- 
equate experimental data are available, first-principles 
calculations provide the necessary data for fitting the 
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TABLE III: Parameters for the Pt-based alloy potential. 



Parameter 


Pt-Pd 


Pt-Cu 


Pt-Ag 


Pt-Au 




1.701 


1.734 


1.820 


1.588 


Dm (eV) 


1.259 


0.770 


0.760 


0.792 


Rm (A) 


2.448 


2.517 


2.555 


2.430 



cross-interaction potentials, such as the formation en- 
thalpy and the lattice constant of the alloys, and, as 
a result, more accurately parameterized potentials can 
be constructed for individual alloys. However, we are 
rarely aware of the development of potentials that were 
optimized to reproduce the properties of various binary 
Pt-based alloy systems. In this study, we optimized the 
adjustable parameters for the cross-interaction based on 
both the formation enthalpy and the lattice constant of 
the alloys. (For cross-interaction, 5 was set to zero.) The 
determined potential parameters for four Pt-based alloys 
are shown Table IIIII 

Figure 13] shows the equilibrium lattice constants of four 
Pt-based alloys as a function of the Pt mole fraction to- 
gether with the published experimental data^"— . (For 
Ag-Pt, we could not find experimental data on the mix- 
ing enthalpy.) Because the Ag-Pt system forms an in- 
termetallic compound, experimental data were not com- 
pared for Ag-Pt. In the case of Pd-Pt, Cu-Pt, and Au- 
Pt, the present model showed reasonable agreement with 
the experimental results. Note that the Cu-Pt system 
showed somewhat positive deviation from Vegard's law, 
whereas the other alloys all showed slightly negative de- 
viation from Vegard's law. 

Figure m shows the heat of formation for Pd-Pt, Cu-Pt, 
Ag-Pt, and Au-Pt bulk alloys. The calculated heats of 
formation data for disordered solid solutions were com- 
pared with the experimental data. Our present model 
showed good agreement with the experimental data and 
the empirical model predictions. Interestingly, the Pd- 
Pt and Cu-Pt systems showed negative heat of formation 
behavior, whereas the Au-Pt and Ag-Pt system showed 
positive heat of solution behavior in the bulk phase. 
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FIG. 3: (a) Lattice parameters and (b) their variation of the 
departure from Vegard's law for the Pd-Pt, Cu-Pt, Ag-Pt, 
and Au-Pt bulk systems as a function of the Pt concentration. 
The lines are the results of the present EAM model, and the 
solid symbols represent the experimental data. 



III. COMPUTATIONAL METHOD 

A. Initial Structures of the Nanoparticles 

It is well known that the metallic nanoparticles present 
competitive structural motifs, such as the icosahedron, 
decahedron, cuboctahedron, and truncated octahedron 
(TOh)^!^. In our atomistic simulations, a truncated 
octahedron (TOh) nanoparticle of 1654 atoms was used 
as an initial atomic configuration, which corresponds to 
~3 to 4 nm in size, depending on its lattice distance. 
The TOh shape of the nanoparticles was chosen because 
it is the experimentally observed shape for this size of fee 
metal nanoparticle^. As shown in Fig.[5l the TOh shape 



is closed-shell structures having 6 {100} facets, 8 {111} 
facets, 12 {111}/{111} edges, 24 {111}/{100} edges, and 
24 vertices and the outermost surface layer is composed 
of 582 atoms (35%). Although we investigated the struc- 
tures of various mole fractions, the composition of 50 at% 
Pt was chosen as an example condition to simplify the 
discussion. With the composition fixed at 1:1 stoichiom- 
etry, it appears that any combinations of the core-shell 
structure are possible. A starting compositional configu- 
ration is typically generated in a random way according 
to the nominal composition of the bulk alloy [Fig. [5] (b), 
although we also used alloy nanoparticles with some ini- 
tially ordered configuration for comparison [Fig.[S](c) and 
(d)]. All of the atoms were fully relaxed with molecular 
statics. 
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FIG. 4: Heat of formation for tiie Pd-Pt, Cu-Pt, Ag-Pt, and 
Au-Pt bulk systems as a function of Pt concentration. The 
solid lines are the results of the present EAM model for a dis- 
ordered soUd solution (random mixing). The solid symbols 
represent the experimental measurement of the heat of for- 
mation (Refs. [SJ and l3ll) and the dashed lines represent the 
model equations that are derived from the experimental data 
and the phase diagrams (Refs. and [3^. 




FIG. 5: Initial atomic configurations: (a) a monometallic 
nanoparticle of fee structure and bimetallic nanoparticles of 
(b) random mixing, (c) LIq and (d) Lli ordered structure. 



B. Monte Carlo Simulations of the Nanoparticles 

As a first step to investigate the structural behavior 
of alloy nanoparticles, we focused on the energetics. For 
nanoparticles, the prediction of energetically stable struc- 
tures is regarded as an optimization problem of high com- 
plexity because it is a problem of searching for the lowest 
energy configuration in the solution space composed of all 
feasible atom arrangements. Although it is not perfect 
for global optimization, Monte Carlo (MC) simulation 
is the most popular method for this purpose^. In this 
study, we performed a series of Monte Carlo simulations 
based on the Metropolis algorithm, where not only the 
number of total atoms but also that of each element in the 
alloy is kept constant: The zero-temperature minimum- 
energy configurations of alloy nanoparticles were pre- 
dicted by optimizing their total energies with respect 
to their compositional exchange. Because the prediction 
quality of a Monte Carlo simulation might depend on 
the initial condition, we adopted various initial configu- 



rations, including a random mixture and some ordered 
structures, as shown in Fig. [5l The zero-temperature 
Monte Carlo simulations were performed for at least 10^ 
Monte Carlo steps (MCS) for the optimizations, where 
one MCS corresponds to N attempted exchanges and 
N is the number of atoms in the nanoparticle. For the 
simulations, we used a hybrid scheme where structural 
relaxation using the conjugate-gradient method was im- 
plemented in the trials of compositional exchange^l. To 
avoid being trapped in a local minimum, we also per- 
formed a series of simulated annealing (SA) optimizations 
combined with the Monte Carlo simulations-^. Here, the 
simulated annealing procedure was performed from T = 
1000 K to K at a cooling rate of of - lO^^ K/MCS. 



IV. RESULTS AND DISCUSSION 

A. Structure of TOh Alloy Nanoparticles 

Although numerical variation was present in the fi- 
nal binding energy of the optimized structures, both the 
simulated annealing and the simple Monte Carlo simu- 
lations produced nearly identical views of the bimetallic 
nanoparticle structures. Also, the optimized structures 
were relatively insensitive to the initial structures, prob- 
ably because of the low transition barrier for composi- 
tional configuration change from a random structure to 
the optimized structure. Figure [6] shows the lowest en- 
ergy configurations of the 1654 atoms of TOh nanopar- 
ticles for the four alloy systems. Of course, the bimetal- 
lic nanoparticles show different compositional configu- 
rations depending on the component metals: (1) The 
Pd-Pt nanoparticle exhibits a typical LIq intermetallic 
compound structure with weak surface segregation; (2) 
The Cu-Pt nanoparticle exhibits a onion-like multi-shell 
structure with Cu on the outer surfaces; (3) The Au- 
Pt and Ag-Pt nanoparticles exhibit Au@Pt and Ag@Pt 
core-shell structures as minimum energy states. 

The inner structure can be analyzed more clearly in 
Fig. [3 where the composition profiles of Pt in each of 
the shells of the nanoparticles were plotted. In the Pd- 
Pt nanoparticle, Pd and Pt atoms are comparably dis- 
tributed over the shells, indicating ordered structures. 
The Cu-Pt nanoparticle exhibits an alternating distribu- 
tion of Cu and Pt atoms over the shells: The 1st, 3rd, 
and 5th shells are Cu-rich, whereas the 2nd, 4th, and 6th 
shells are exclusively occupied by Pt atoms. Note that 
after peeling off the Cu-rich outer shell, we can obtain a 
perfect Pt surface. In Au-Pt and Ag-Pt nanoparticles, Pt 
atoms exist only on the core region and are completely 
surrounded by the An or Ag shell (Au@Pt and Ag@Pt 
core-shell nanoparticles) . 
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cross-sectional drawing surface (100) facet tilt view 

FIG. 6: Predicted atomic structures of (a) Pd-Pt, (b) Cu-Pt, 
(c) Au-Pt, and (d) Ag-Pt nanoparticles. The first, middle, 
and last columns of figures show the cross-section, surface, 
and tilt view of the nanoparticles, respectively. The TOh 
nanoparticles are composed of 1654 atoms (M827Pt827, where 
M stands for Ag, Au, Cu, and Pd). Pt atoms are in grey 
(light), Pd in black, Cu in blue, Au in green, and Ag in red 
(dark). 
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FIG. 7: Pt compositional profile in the shells of the 1654 
atom TOh alloy nanoparticles. The inset shows seven shells 
of the nanoparticle in different colors: the first shell number 
denotes the surface, the second shell denotes the subsurface, 
and so on. 



B. Equilibrium Structure of Bulk Alloy Phases 

Some material properties that seem to be deeply re- 
lated to alloy structures are already shown in Fig. [T] (sur- 
face energy). Fig. [2] (cohesive energy and melting point), 
and Fig. |4] (heat of random mixing) in Set. HIl In ad- 
dition, the calculated heats of formation for three dif- 



ferent structures of MgoPtso (M 50% - Pt 50%, where 
M stands for Ag, Au, Cu, and Pd) bulk alloys, includ- 
ing one disordered and two ordered phases, are given in 
Fig.m As previously shown in Fig.|31 Pd and Cu exhibit 
negative heats of mixing with Pt, whereas Au and Ag 
exhibit positive heats of mixing behavior with Pt in the 
bulk phases. At the particular composition of 1:1, the 
tendency for ordering is strong in both Pd-Pt and Cu-Pt 
alloys2^. For example, our potential model predicts that 
the Pd-Pt bulk alloy clearly favors the intermetallic com- 
pound phase of the LIq ordered structure as a minimum 
energy state. The Cu-Pt alloy shows even larger nega- 
tive heat of mixing and CusoPtso favors the intermetallic 
compound phase of the Lli ordered structure over ran- 
dom mixing or the LIq structure. Interestingly, CuPt is 
known as the only metallic alloy system that forms in 
the Lli structure'*^. This structure consists of alternat- 
ing fee (111) layers of Cu and Pt, in contrast with the LIq 
structure of alternating (001) planes of atoms, as shown 
in Figs.[S](b) and (c). 

Our EAM alloy model predicts that the orders of for- 
mation energy for the three bulk alloy structures are 
A7?(Lli) > Ai/ (disordered) > Ai?(Llo) for Pd-Pt and 
AT? (disordered) > AiT(Llo) > Ai/(Lli) for Cu-Pt, re- 
spectively. Note that the predicted lowest energy crystal 
structures for PdsoPtso and CusoPtso bulk alloys are in 
excellent agreement with the experiment c'^^ ' ' — '. In con- 
trast, the formation energy of mixed structures is less 
important for Au-Pt and Ag-Pt, because none of them is 
an equilibrium structure over the separated phase. Re- 
gardless, because the bulk structure is deeply related to 
the inner structure of the nanoparticles, accurate repro- 
duction of the bulk alloy structures is critical for reason- 
able prediction of the optimized structures of the Pd-Pt 
and Cu-Pt nanoparticles. 

C. Structural features of Pd-Pt and Cu-Pt 
nanoparticles 

As depicted in Fig. [HI the tendency for ordering is 
strong in both the Pd-Pt and Cu-Pt alloy nanoparticles. 
The Pd-Pt nanoparticle forms a mixed alloy of ordered 
structure, as does the bulk PdPt alloy. However, the 
Cu-Pt nanoparticle forms a unique crystal structure that 
cannot be observed in bulk structures: it forms an onion- 
like structure (with surface segregation of Cu) rather 
than the LIq or Lli structures. To analyze the bond 
characteristics of the alloys, we analyzed the structural 
properties with the Warren-Cowley chemical short-range 
order (CSRO) parameter defined as^ 

where Nab is the nearest coordination number of B 
atoms around an A atom, N is the total coordination 
number in the nearest-neighbor shell, and Cb is the 
atomic concentration of B. 
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FIG. 8: Three types of binary alloy structures: (a) disordered 
solid solution, (b) Llo, and (c) Lli crystal lattice, (d) Heats 
of formation for the three alloy crystal types in the Pd-Pt, 
Cu-Pt, Au-Pt, and Ag-Pt alloy systems. 




PdPt CuPt AuPt AgPt 

(^) Alloy nanoparticle 

FIG. 9: Local atomic configurations in (a) Llo and (b) Lli 
crystal lattices; (c) Warren-Cowley chemical short-range or- 
der (CSRO) parameters for the four alloy nanoparticles. 



This parameter is a useful measure of the chemical 
affinity and represents the degree of tendency for order- 
ing or clustering: CSRO values range between -1 to 1 and 
positive CSRO parameter values indicate clustering or 
phase separation and negative values indicate strong A- 
B bond and chemical ordering. As depicted in Figs. EUb) 
and (c), the LIq structure consists of alternating fee (001) 
planes of Cu and Pt, whereas the Lli structure consists 
of alternating fee (111) planes of atoms. These atomic 
arrangements of different alloy lattice types can also be 
distinguished by the CSRO values. For example, atoms 
in the LIq crystal structure have 4 like-atom neighbors in 
the same (001) plane and 8 (=4-1-4) unlike-atom neigh- 
bors in two nearby (001) planes as shown in Fig. |9l[a). 
So the local CSRO value should be 1/3 for the LIq crys- 
tal structure. Likewise, atoms in the Lli crystal struc- 
ture have 6 like- atom neighbors in the same (111) plane 
and 6 (=3-1-3) unlike-atom neighbors in the nearby (111) 
planes, as shown in Fig. |9jb), and the local CSRO value 
goes to zero, regardless of their element type. 

Figure M,c) shows the CSRO parameter for the four 
alloy nanoparticles presented in Fig. [5] Whereas the 
Au-Pt and Ag-Pt nanoparticles exhibit positive CSRO 
parameters with values close to 1, indicating strong sep- 
aration, the CSRO of the Pd-Pt nanoparticle is approx- 
imately —0.33, which is nearly the ideal value of the 
bulk Llo structure, and the inset to the figure indicates 
nearly identical ordering inside the nanoparticle for the 
Llo structure. However, the CSRO of the Cu-Pt onion- 
like nanoparticle is very close to zero, which is the ideal 



value for bulk Lli structures. Although the apparent 
atomic configurations appear to be quite different from 
each other, the onion-like structure is, however, very sim- 
ilar to Lli in terms of the chemical short-range-order 
(CSRO). Because they are both ordered structures of 
multi-layers or multi-shells, the overall chemical bond 
status is expected to be very similar. 

The transition of the lowest energy structure of certain 
alloy nanoparticles from the Lli structure to an onion- 
like structure is very interesting, and we expect that it is 
probably due to surface effects. To understand the role of 
surface energy in nanoparticle structures, we summarized 
the (001) and (111) surface energies of various relevant 
Pd-Pt and Cu-Pt alloy structures that can appear in TOh 
nanoparticles, as shown in Table IIVI All of the surface 
energies were calculated based on our present EAM po- 
tentials using a slab of more than 100 atomic layers under 
a 2-dimensional periodic boundary condition. 

For Pd-Pt, both the (100) and (111) surface energies 
of Pd are comparable to those of Pt. The surface en- 
ergies of the alloy Llo phase, which is the equilibrium 
state in bulk, are also comparable to those of pure Pt or 
Pd (within less than a 10% difference). The negligible 
dependence of surface energy on alloy composition may 
be attributed to the similar atomic sizes and bonding 
characteristics of Pd and Pt, such as cohesive energy and 
melting point. As shown in Fig. [SJ ideal TOh nanopar- 
ticles of Llo structure can have several different types of 
facets, including (100) facets of pure Pt or Pd and (111) 
facets of alloy composition. Because of the comparable 
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TABLE IV: Comparison of the surface energies for various 
surface types. The numbers in parentheses are relative values 
with respect to the surface energy of pure Pt. 



Bulk structure 


Surface enerj 
(001) 


;y [mJ/m^] 
(111) 


Pt (fee) 


1891 (1.00) 


1580 (1.00) 


Pd (fee) 


1815 (0.96) 


1602 (1.01) 


Cu (fee) 


1435 (0.76) 


1266 (0.80) 


PdPt (Llo) 


1843 (0.97) 


1631 (1.03) 


CuPt (Lli) 


1846 (0.98) 


1967 (1.25) 



surface energy properties between those facets, the Pd- 
Pt nanoparticles can maintain a similar structure to the 
bulk Llg phase with only a small surface modification, 
even at this nanoscale. 

For the Cu-Pt system, the (100) and (111) surface en- 
ergies of Cu are much lower than those of Pt: the surface 
energy ratio of Cu with respect to that of Pt is approx- 
imately 0.76 for the (100) surface and 0.80 for the (111) 
surface. However, the surface energies of the alloy Lli 
phase, which is the equilibrium state in bulk, are compa- 
rable or even higher than those of Pt, let alone Cu. In 
particular, the (111) surface energy of the CuPt alloy Lli 
phase, one of the dominant facet types in the TOh shape, 
is higher than that of Cu by ^ 50%. The remarkable in- 
crease in surface energy when forming an alloy structure 
may be partially attributed to the large atomic size mis- 
match between Pt and Cu (by approximately 10%) that 
can produce surface waves. Because the surface energy 
tends to follow the bulk binding energy, the large nega- 
tive heat of formation may also contribute to the large 
surface energy of the alloy composition. In any case, 
the surface energy difference between the two elements, 
which is not critical in bulk mixtures, becomes a crucial 
factor in nanoparticle formation: The more than 20% 
difference in surface energy between Cu and Pt/Cu-Pt 
leads to the pure Cu (100) and Cu-rich (111) facets, re- 
sulting in the onion-like structures, as shown in Fig. [51 
This mechanism can operate effectively in minimizing the 
surface energy, especially when one of the component 
species can take advantage of a significantly lower sur- 
face energy. Due to the interplay of chemical ordering 
and the surface energy effect, we now believe that an al- 
loy of Lli ordered structure in bulk tends to favor an 
onion-like multi-shell structure when it forms a nanopar- 
ticle. To this end, some alloy systems have been reported 
to exhibit onion-like structures in previous reports of MD 
simulationai^ii^. However, we found that the structure 
of the alloy nanoparticle is very sensitive to the empir- 
ical potentials adopted in the simulation and care must 
be taken in the interpretation. 



D. Core-Shell Structures (AuOPt and AgOPt) of 
Au-Pt and Ag-Pt nanoparticles 

Although Pt, Au, and Ag are fee metals, and the mis- 
match of the lattice constants between Pt and Au or Pt 
and Ag is relatively small (~ 4%, see Table iHl) . the Au- 
Pt and Ag-Pt alloy systems are thermodynamically im- 
miscible, exhibiting large positive heat of mixing values. 
Referring to the phase diagrams of bulk alloy systems, 
we know that Ag and Pt form a peritectic system and 
that Ag and Pt form a miscibility gap. Also, the phase 
diagrams indicate that the mutual solubility of the Au- 
Pt and Ag-Pt system is almost zero or less than 4% at 
room temperature^fl. Therefore, a bulk solid alloy of 1:1 
composition would decompose into the fee fee ai and a2 
phases of limited solid solubility. Regarding the forma- 
tion of nanoparticles, our simulations predict that the 
Au-Pt and Ag-Pt alloy nanoparticles energetically favor 
Pt-core-Ag-shell and Pt-core-Au-shell atomic structures. 
In addition to the phase separation tendency due to the 
positive heat of mixing, the formation of the core-shell 
structure can be interpreted by the surface energy effect. 
As shown in Fig. [U the surface energies of both Au and 
Ag are only approximately 50 to ~ 60% of that of Pt. 
Because of the differences in surface energy between Pt 
and Au/Ag, Pt always forms the core, whereas Au or Ag 
forms an outer shell to minimize the surface energy. Our 
prediction of core-shell structures (Au@Pt and Ag@Pt) 
for Au-Pt and Ag-Pt nanoparticles agrees well with the 
first-principle calculation"'^ and Monte Carlo studies^—, 
although the size and shape conditions are slightly dif- 
ferent. 



V. CONCLUSION 

For systematic design of bimetallic nanoparticle cata- 
lysts, we developed a unified EAM potential model for 
various Pt-based binary alloys, such as Pd-Pt, Cu-Pt, 
Au-Pt, and Ag-Pt. After verifying that our alloy poten- 
tials are reliable enough to accurately reproduce the sur- 
face energy and heat of formation of alloy mixtures, we 
performed a series of simulated annealing Monte Carlo 
optimizations on TOh-shaped bimetallic nanoparticles 
consisting of 1654 atoms. Our simulations quantify the 
energetically favorable atomic arrangements of the four 
Pt-based alloy nanoparticles: the Lip intermetallic com- 
pound structure for the Pd-Pt system, an onion-like 
multi-shell structure for the Cu-Pt system, and core-shell 
structures (Au@Pt and Ag@Pt) for the Au-Pt and Ag-Pt 
systems. The main features of these optimized nanopar- 
ticle structures were analyzed on the basis of the rela- 
tionship with their bulk alloy structures. The equilib- 
rium structure of bimetallic nanoparticles could be in- 
terpreted by the interplay of material properties, such as 
the surface energy and heat of formation. We believe that 
understanding the energetically favorable structures can 
be a first step to the new design of nanoparticle struc- 
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tures. To this end, the relationship between the material 
properties and the resulting nanoparticle structures is of 
primary importance; other parameters, such as the size, 
composition, and temperature effect will be investigated 
in future studies. 
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